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Abstract 

The Gumowski-Mira map is a 2-dimensional recurrence relation that provide a large 
variety of phase space plots resembling fractal patterns of nature. We investigate 
the nature of the dynamical states that produce these patterns and find that they 
correspond to Type I intermittency near periodic cycles. By coupling two GM maps, 
such patterns can be wiped out to give synchronised periodic states of lower order. 
The efficiency of the coupling scheme is established by analysing the error function 
dynamics. 
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1 Introduction 



Non linear recurrence relations model many real world systems and help in 
analysing their possible asymptotic behaviour as the parameters are varied [1]. 
Here we analyse one such recurrence map defined through a 2-dimensional 
difference equation called Gumowski-Mira transformation [2]. This system in- 
volves a non linear function f{x) that is advanced in time by one iteration, 
in the y equation. It has been reported that computer simulations of this it- 
erative scheme [3] give rise to a variety of 2-dimensional images in the phase 
plane (x, y) called GM patterns that resemble natural objects like star fish, 
jelly fish, wings of a butterfly, sections of fruits, flowers of varying shapes etc. 
Apart from a mere fascination in creating these patterns, one can implement a 
scheme like Iterative function scheme for producing 2-dimensional fractals or 
fractal objects of nature using this map for different values of the parameters 
involved. Moreover, the final pattern is found to depend very sensitively on the 
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parameters, a feature which can be exploited in decision making algorithms 
and control techniques in computing and communication. 

One motivation in this work is to analyse the underlying dynamics behind 
the formation of these patterns and their possible bifurcation sequences. In 

this context, we would like to mention that the period doubling sequences 
and chaos doubling phenomenon in this system has been reported earlier for 
a chosen set of parameter values [4]. We are concentrating on another set 
that was used in [3] for generating GM patterns. Here the dynamics involved 
is different in the sense that the intermittency route and odd period cycles 
predominate and the route to chaos is via formation of quasiperiodic bands 
and band mergings. The detailed evolution of the basin structure during the 
scenario until onset of chaos and final escape forms part of our work. The 
scaling behaviour near intermittency, stability of the lowest periodic cycles are 
also analysed. Moreover our work provides evidence for periodic window and 
repeated substructures in the bifurcation scenario which show self similarity 
with respect to scaling in the parameter. This can therefore account for the 
sensitive dependence of the asymptotic dynamical states or phase portraits on 
the relevant parameter of the system. 

It has been reported that [5] the lowest order periodic orbits of the GM map 
are non parametric in nature and hence a straight forward application of OGY 
scheme to control chaos in the system is not desirable. Hence this is achieved 
through a proxy system in [5]. We introduce a simple and straight forward 
linear and mutual coupling scheme to achieve the same result. Moreover by 
coupling two such maps synchronised states with low periodic cycles can be 
produced. Periodicity with synchronisity is found to be quite robust in the 
system even against considerable perturbations. A detailed analysis of the 
error function dynamics is carried out to estimate the average synchronisation 
time (ti) and the stabilisation time {T2) after a perturbation is applied. 

In this paper. Section 2 introduces the GM map and its possible stable dynam- 
ical states. Section 3 gives the results of a detailed numerical analysis of the 
predominant bifurcation sequences in the map. The intermittency and other 
routes to chaos are given in Section 4. In Section 5, the control of chaos and 
synchronisation of the coupled maps is discussed. The discussion of the main 
results are given in Section 6. 



2 Dynamical States of GM map. 

The transformation of Gumowski-Mira map that has been studied in computer 
simulations can be written as a 2-dimensional recurrence relation defined in 
the X — y plane as 
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Xn+1 = ?/n + a(l - hyl)yn + f{Xn) 

Vn+l = ~Xn + f{Xn+l) 



(1) 



with ,f{xn) = fJ'Xn + ^|j^~2^ x^. It should be noted that the function /(x) in 
the y equation is advanced in time by one iteration which provides richer 
and complex dynamical behaviour. The nature of variation of the non linear 
function /(xq) as well as f{xi) for different // values in the range —0.3 to 
+0.0039 with a = 0.008 and b = 0.05 is shown in figs la and lb. Here Xq is the 
starting seed value and Xi is the first iterate. These figures give an indication 
of how different the functions in the x and y equations arc and how fast and 
differently they change as iterations proceed. As any stable periodic cycle is 
reached asymptotically, these functions are found to level off. 



The one cycle fixed points of the system is determined by solving the equations 



X ^ y + a{l — hy'^)y + + 



2(l-;u)a;2 
l + a;2 



y^{n-l)x + 



2(1 - n)x^ 
1 + 



(2) 



using MAXIMA and the only real and bounded solutions of interest are (0, 0) 
and (1, 0), which are unstable in the relevant parameter regions of our study. 
The most prominent elementary cycle here is a 4 cycle born by saddle node 
bifurcation. The stability of this 4 cycle can be established numerically by 
calculating the cigcn values of the Jacobian M for chosen values of /i and 
the corresponding elements of the 4 cycle. Thus if {Sxn, Syn) are the small 
perturbations to the cycle elements, we have 



^SXn+l 



M 



(3) 
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where 

Mi2 =1 - 2a6?/' + a(l - 6?/') 

4a;3(l-/x) 4a;(l-//) 

_4xf(l-M) , 4x(l-M) , ,A Z',, , _ L„2N , 2x^(1-/.) , ^, \3 



M2i^-i + //|- ,/ 7 + ;_r +/^l- 



(! + (?/ + a|/(l - + + ^^)2)^ 

4(1-/") + ^ + (y + ay(l - 6^^) + + x/.) 

1 + + ay(l - by^) + + 



(4) 



M22 =(1 - 2a6?/2 + _ 6^2))^- 

4(1 - 2aby' + a(l - by')){l - fx){y + ayjl - by') + '-^^ + xj^f ^ 

(1 + (y + ay{l - by^) + + xf^rf 

4(1 - 2aby' + a(l - 6y^))(l - fx){y + ayjl - by') + '-^^ + x//) 
l + (z/ + ay(l-6?/2) + ?^JM + ^^)2 

The eigen values of M evaluated for a few typical jj, values with a = 0.008; 
b = 0.05 in the 4 cycle window are given in 1 . 

Table 1 

The Eigen values of the stability matrix for different ^ values in the 4 cycle window 
|£^| < 1 in all the above cases indicating stability of the concerned 4 cycle [6] 





El 


E2 


-0.0964 


0.6384 


-0.6346 


-0.095 


0.6451 


-0.6415 


-0.005 


0.9371 


-0.9371 


0.0002 


0.8169 


-0.8169 


0.0003 


0.7841 


-0.7841 



3 Bifurcation Sequence. 



We have mentioned that the asymptotic dynamical states of the GM map 
depends sensitively on the parameter /x. In this section, we numerically inves- 
tigate in detail, the possible states and their bifurcation patterns as is varied. 
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We find that the bounded interval for the map lies in the interval [—1,1] for 
a = 0.008 and b = 0.05. Moreover the most prominent elementary cycle of 
periodic behaviour is 4 which occurs in many intermittent windows of /j, in 
the bifurcation diagram. The full scenario is given in Fig 2a which is mostly 
dominated by broad windows of odd cycles hke 7, 11 etc. Here out of 10,000 
iterate 9000 are discarded as initial transients and the next 1000 are plotted. 
The specific regions of the windows of such cycles are zoomed and reproduced 
in Fig 2b, 2c and 2d. 

In general these windows of periodic cycles born by tangent bifurcations or 
saddle node bifurcations at their left ends exhibit intcrmittcncy behaviour in 
their iterates. The transition to chaos takes place for small increase of values 
of when the periodic cycle becomes unstable giving rise to quasiperiodic 
bands which become chaotic and merge together. We illustrate this for the 
7 cycle window. Near the left end of a periodic 7 cycle window Fig 3a gives 
the — n plot with ^ = —0.2734 for 6000 iterations showing intermittently 
laminar and chaotic behaviour. Fig 3b gives the corresponding x — y plot. It 
is interesting to note that the periodic window structure shows self similarity 
in their substructures with respect to scaling in /i. This is evident from the 
plots in Fig 3c and Fig 3d where a 7 cycle and a 22 cycle patterns repeat at 
two levels of zooming in the parameter range. This recurring self similarity 
with the consequent intermittency in the beginning of each periodic cycle and 
quasi-periodicity at its end, accounts for the variety and richness of the GM 
patterns in the x — y plane with its highly sensitive dependence on fj,. 



4 Intermittency &c Transition to Chaos. 

The nature of intermittency before the birth of a periodic cycle is analysed by 
calculating the average life time of the laminar region (/) in the a;„— n plots. For 
one such region near /i — —0.312501 before a ten cycle window the variation 
of (l) as a function of Ifi — fid is studied (Fig 4). /x^ is the critical value when it 
becomes chaotic. We can write (/) ~ — yud*^. v then defines the scaling index 
that helps to identify the type of intermittency. Here iic = —0.312498 and 
u — —0.49 so that the intermittency is Type I [7] which is usually associated 
with saddle node bifurcation. 

For further work we concentrate mainly on the stable window of the 4 cycle 
in the range -0.096 < /i < 0.0003695. 

The phase space during intermittency before the stabilisation of the 4 cycle, 
the 4 cycle region, the quasi periodic band region with 4 bands and final 
merging of bands to chaos are shown in Fig 5a-d. It is clear that the transition 
to chaos is via quasi-periodicity to chaotic bands and band merging. The 
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precise transition point is in this case ji — 0.000369. 

The evolution of the basin structure in the phase space during these transitions 
is also studied in detail and shown in Fig 6a-c. 

The last figure 6d gives the basin boundary between bounded and escape 
regions, where escape takes place beyond ji — 



5 Synchronisation and Control of Chaos. 

We try a linear mutual coupling of two GM maps as a control mechanism, by 
which chaotic and intermittency regions can be targeted to stable low periodic 
regions. 

The corresponding equations are 

Xln+l = yln + a{l - byll)yln + fixln) + Sx2n 
y^n+l = -Xln + f{xln+l) 

x2n+i = y2n + a(l - by2l)y2n + /(x2„) + exln 

y2n+i = -x2n + f{x2n+i) (5) 

where e is the coupling parameter. Wc analyse this system numerically for 
fi = —0.31 and a = 0.008, b = 0.05, both the maps individually exhibit chaotic 
behaviour. After coupling with 6 = 1 and initial values xIq = 0.1, ylo = 
for the first map and x2o = 0.2 and y2o = 0.1 for the second map, both the 
maps are found to settle to 5 cycles. Here control of chaos and periodicity are 
achieved eventhough synchronisation is absent. 

We are also able to synchronise two chaotic maps to lower periodicities us- 
ing this coupling scheme. Eventhough much work is reported in the area of 
synchronisation [8] , majority of these works are in continuous systems. How- 
ever there are a few specific cases reported in the context of two dimensional 
map [9,10]. Here we achieve synchronisation in two coupled discrete maps of 
the GM type for /j, = —0.39 with e = 0.7. Total synchronisation is seen in both 
X and y as the coupled system settles to identical 4 cycles. This is shown in 
Fig 7a. 

The analysis is continued by developing the error function dynamics by defin- 
ing the error in the x and y values of the two synchronisating systems as 
= {xl — x2) and = (yl — y2). Then their dynamics develops through the 
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following set of maps. 



/' ■ 



2{1 - ^){xln + X2r. 



+ 



1 + Xl2)(l + X22) 

1 + a - ab{yll + ?/l„|/2 
2(1 - ^)(xl„+i 



(1 + + X2^_j_;^) 



(6) 



The set of equations (5) and (6) are evolved together until — >• and 
[11], when total synchronisation is said to be achieved. 

From this the time for synchronisation is computed for 10 sets of initial values 
and the average is found to be (ri) = 502 for the typical values mentioned 
above. 

The robustness of the mechanism is tested by applying a perturbation and 
the stabilisation time is calculated as the time to reach synchronisation again 
after it is disturbed. This is also repeated for 10 initial values and the average 
time {T2) is found to be 473. 

Fig 7b is the — n plot showing ri and T2 values for a perturbation of Xn = 
Xn + 10 given after 2000 iterations. A similar plot is obtained for the error 
function in y also. 



6 Conclusion 



In this work we analyse the dynamics of the large variety of interesting and 
hvely patterns exhibited by the GM map and find that they are attr actors of 
the system in the phase plane in the neighbourhood of higher order cycles. 
The temporal behaviour then corresponds to intermittency of Type I with an 
exponent ~ —0.5. 

The dependence of the phase space structure of the patterns on minute changes 
in the parameter makes it a useful tool in decision making algorithms. This 
sensitive dependence is due to the recurring periodic and self similar substruc- 
tures in the bifurcation scenario, each with its own intermittency, periodicity, 
quasi periodic band and merging of bands leading to chaos. The nature and 
characterisation of these patterns and occurrence of crises related phenomena 
near them are being studied and will be reported else where. 

The coupled states of these maps are either synchronised or they stabilise to 
attractors of the same periodicity, without amplitude synchronisation. A de- 
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tailed analysis of the error function dynamics is carried out to estimate average 
time for synchronisation and stabilisation time after applying a perturbation. 

The relative efficiency of different schemes of coupling in synchronising 2 such 
maps to chosen dynamical states is intended as further study. 
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Fig. 1. Variation of a) the function f{x) b) its first iterate for three different values 
of II -0.3, -0.16 and 0.0039 
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-0.2 -O.IS -0.16 -0.L4 -0.12 -0.1 -O.OS -0.06 -O.l 



O.QL -O.OOS -0.006 -0.004 -0.002 



(c) -0.2 <n< -0.0505 (d) -0.01 < /x < -0.00095 

Fig. 2. a) Bifurcation scenario in the range jj, = —0.2 to +1, for a = 0.008, b = 0.05. 
The windows of periodic cycles are zoomed for details in b, c and d as indicated. 
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(c) 



(d) 



Fig. 3. The intermittency behaviour before the stabilisation of a 7 cycle is shown in 
the Xn-n plot in (a) and corresponding phase portrait in (b) for jjl = —0.2734. Note 
that this is the typical behaviour near the onset of each periodic window. The self 
similar and repeating substructures inside this stability window is shown in Figs 3c 
and 3d where the 7 and 22 cycles are seen to recurr. (The rectangle shown in 3c is 
zoomed in 3d). 
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(a) (b) 



Fig. 4. Type I intermittency near the onset of a 10 cycle a) the laminar and irregular 
behaviour for fi = —0.312501 b) the scaling of average laminar region (/) as function 
of l/Lt — /Xc| where = —0.312498. The scaling index in this case is —0.49. 



(a) (b) 




(c) 



(d) 



Fig. 5. The phase portraits for jj, values in the region of the 4 cycle window a) 
intermittency region (p = —0.1199) b) stable 4 cycle {p = —0.005) c) quasiperiodic 
4 lands (/x = 0.0003695) and d) chaotic attractor {/j, = 0.001) 
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(c) (d) 

Fig. 6. The basin of attraction for the different dynamical states inside the 4 cycle 
window. Here black region corresponds to the basins of the stable 4 cycle as the 
/J. value increases from intermittency towards periodic cycle, a) ji = —0.0815 b) 
II = —0.081 c) /X = —0.08. In (d), the black region corresponds to bounded chaos 
while white region that of escape for /x = 1.01. 
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a I 1 1 < < 1 1 1 

-8-6-4-2 2 4 6 

(a) The synchronised 4 cycle of two 
coupled GM maps with linear and 
mutual coupling for coupling param- 
eter £ = 0.7, in eqn. (5). 

Fig. 



-30 I < 1 1 < < 1 1 < 1 1 

a 1000 2000 3000 4000 5000 EOOO 7000 SOOO 3000 LOOOO 
n 

(b) The synchronisation time ri and 
the stabilisation time after perturba- 
tion T2 are shown in the — n plot 
of the coupled system 
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